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We introduce the concepts of participant triangularity and triangular flow in heavy-ion collisions, 
analogous to the definitions of participant eccentricity and elliptic flow. The participant triangularity 
characterizes the triangular anisotropy of the initial nuclear overlap geometry and arises from event- 
by-event fluctuations in the participant-nucleon collision points. In studies using a multi-phase 
transport model (AMPT), a triangular flow signal is observed that is proportional to the participant 
triangularity and corresponds to a large third Fourier coefficient in two-particle azimuthal correlation 
functions. Using two-particle azimuthal correlations at large pseudorapidity separations measured 
by the PHOBOS and STAR experiments, we show that this Fourier component is also present 
in data. Ratios of the second and third Fourier coefficients in data exhibit similar trends as a 
function of centrality and transverse momentum as in AMPT calculations. These findings suggest 
a significant contribution of triangular flow to the ridge and broad away-side features observed in 
data. Triangular flow provides a new handle on the initial collision geometry and collective expansion 
dynamics in heavy-ion collisions. 



I. INTRODUCTION 

Studies of two-particle azimuthal correlations have 
become a key tool in characterizing the evolution 
of the strongly interacting medium formed in ultra- 
relativistic nucleus-nucleus collisions. Traditionally, the 
observed two-particle azimuthal correlation structures 
are thought to arise from two distinct contributions. 
The dominant one is the "elliptic flow" term, related 
to anisotropic hydrodynamic expansion of the medium 
from an anisotropic initial state [IHS]. In addition, one 
observes so-called "non-flow" contributions from, e.g., 
resonances and jets, which may be modified by their 
interactions with the medium [T0TH2] . 

The strength of anisotropic flow is usually quan- 
tified with a Fourier decomposition of the azimuthal 
distribution of observed particles relative to the reac- 
tion plane [13] . The experimental observable related 
to elliptic flow is the second Fourier coefficient, "i>2-" 
The elliptic flow signal has been studied extensively in 
Au+Au collisions at RHIC as a function of pseudora- 
pidity, centrality, transverse momentum, particle species 
and center of mass energy [3H7J. The centrality and 
transverse momentum dependence of vi has been found 
to be well described by hydrodynamic calculations, which 
for a given equation of state, can be used to relate a 
given initial energy density distribution to final momen- 
tum distribution of produced particles [13]. In these 
calculations, the t>2 signal is found to be proportional to 
the eccentricity, £27 of the initial collision region defined 
by the overlap of the colliding nuclei [15] ■ Detailed 
comparisons of the observed elliptic flow effects with 
hydrodynamic calculations have led to the conclusion 
that a new state of strongly interacting matter with very 
low shear viscosity, compared to its entropy density, has 
been created in these collisions [H] fTBHTH] . 

Measurements of non-flow correlations in heavy-ion 
collisions, in comparison to corresponding studies in p+p 



collisions, provide information on particle production 
mechanisms [19] and parton-medium interactions (101 - 
112] . Different methods have been developed to account 
for the contribution of elliptic flow to two-particle corre- 
lations in these studies of the underlying non-flow corre- 
lations [XQl [19] [22] . The most commonly used approach 
is the zero yield at minimum method (ZYAM), where 
one assumes that the associated particle yield correlated 
with the trigger particle is zero at the minimum as a 
function of A0 after elliptic flow contribution is taken 
out [2Tj. The ZYAM approach has yielded rich correla- 
tion structures at A</> ss 0° and A(j> ~ 120° for different 
Pt ranges [23Tf26] . These structures, which are not 
observed in p+p collisions at the same collision energy, 
have been referred to as the "ridge" and "broad away- 
side" or "shoulder". The same correlation structures 
have been found to be present in Pb+Au collisions at 
y/s^ = 17 A GeV at the SPS [27]. Measurements at 
RHIC have shown that these structures extend out to 
large pseudorapidity separations of At] > 2, similar to 
elliptic flow correlations [25] . The ridge and broad away- 
side structures have been extensively studied experimen- 
tally [T2J [25H2H1 [25] and various theoretical models have 
been proposed to understand their origin [291436J . A 
recent review of the theoretical and experimental results 
can be found in [37]. 

In this paper, we propose that the observed ridge and 
broad away-side features in two-particle correlations may 
be due to an average triangular anisotropy in the initial 
collision geometry which is caused by event-by-event 
fluctuations and which leads to a triangular anisotropy 
in azimuthal particle production through the collective 
expansion of the medium. It was shown that, in the NEX- 
SPHERIO hydrodynamic model, ridge and broad away- 
side structures in two particle correlations arise if non- 
smooth initial conditions are introduced [36] . Sorensen 
has suggested that fluctuations of the initial collision 
geometry may lead to higher order Fourier components 
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FIG. 1: Top: azimuthal correlation functions for mid-central (10-20%) Au+Au collisions at ^/s NN = 200 GeV obtained from 
projections of two-dimensional Arj, A(/> correlation measurements by PHOBOS |19l 125] and STAR |41| . The transverse 
momentum and pseudorapidity ranges are indicated on the figures. Errors bars are combined systematic and statistical errors. 
The first three Fourier components are shown in solid lines. Bottom: the residual correlation functions after the first three 
Fourier components are subtracted. 



in the azimuthal correlation function through collective 
effects [35]. An analysis of higher order components 
in the Fourier decomposition of azimuthal particle dis- 
tributions, including the odd terms, was proposed by 
Mishra et al. to probe superhorizon fluctuations in the 
thermalization stage |39| . In this work, we show that 
the second and third Fourier components of two-particle 
correlations may be best studied by treating the com- 
ponents of corresponding initial geometry fluctuations 
on equal footing. To reduce contributions of non-flow 
correlations, which are most prominent in short pseudo- 
rapidity separations, we focus on azimuthal correlations 
at long ranges in pseudorapidity. We show that the ridge 
and broad away-side structures can be well described by 
the first three coefficients of a Fourier expansion of the 
azimuthal correlation function 



dNP a 



dA0 



2tt 



2V„a cos(nA</>) 



(1) 



where the first component, Via 1 , is understood to 
be due to momentum conservation and directed flow 
and the second component V2A is dominated by the 
contribution from elliptic flow. Studies in a multi-phase 
transport model (AMPT) [40] suggest that not only the 
elliptic flow term, V 2 a, but also a large part of the 
correlations measured by the V3A term, arises from the 
hydrodynamic expansion of the medium. 



1 Note the distinction between V„a and v n . See Eqs. |l0| and |ll| 
for details. 



II. FOURIER DECOMPOSITION OF 
AZIMUTHAL CORRELATIONS 

In the existing correlation data, different correlation 
measures such as R(At], A<f>) [Jj5], Nf(Ar], A<fi) [JJJ and 
l/A r trigdA^/dA0(A?7, A<f>) [M] have been used to study 
different sources of particle correlations. The azimuthal 
projection of all of these correlation functions have the 
form 



C(A(f>) = A 



dA(j) 



B, 



(2) 



where the scale factor A and offset B depend on the 
definition of the correlation function as well as the 
pseudorapidity range of the projection [23] . Examples of 
long range azimuthal correlation distributions are shown 
in Fig. [T] for mid-central Au+Au collisions with different 
trigger and associated particle px selections obtained by 
projecting the two-dimensional correlation functions onto 
the A(j) axis at pseudorapidity separations of 1.2 < Ar] < 
1.9 for STAR data [UJ and 2 < A77 < 4 for PHOBOS 
data [Till [22] . The correlation function data used in this 
study are available at [4"2"H4"4l . Also shown in Fig. [I] 
are the first three Fourier components of the azimuthal 
correlations and the residual after these components are 
taken out. The data is found to be very well described 
by the three Fourier components. 



III. PARTICIPANT TRIANGULARITY AND 
TRIANGULAR FLOW 

It is useful to recall that traditional hydrodynamic 
calculations start from a smooth matter distribution 
given by the transverse overlap of two Woods-Saxon 
distributions. In such calculations, elliptic flow is aligned 
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FIG. 2: Distribution of h 
in Jsnn = 200 GeV Au+Au collisions. 



eccentricity, £2, and (b) triangularity, £3, as a function of number of participating nucleons, N pal t, 



with the orientation of the reaction plane defined by the 
impact parameter direction and the beam axis and by 
symmetry, no V3A component arises in the azimuthal 
correlation function. To describe this component in 
terms of hydrodynamic flow requires a revised under- 
standing of the initial collision geometry, taking into 
account fluctuations in the nucleon-nucleon collision 
points from event to event. The possible influence of 
initial geometry fluctuations was used to explain the 
surprisingly large values of elliptic flow measured for 
central Cu+Cu collision, where the average eccentricity 
calculated with respect to the reaction plane angle is 
small [SJ. For a Glauber Monte Carlo event, the minor 
axis of eccentricity of the region defined by nucleon- 



nucleon interaction points does not necessarily point 
along the reaction plane vector, but may be tilted. The 
"participant eccentricity" [3 [45] calculated with respect 
to this tilted axis is found to be finite even for most 
central events and significantly larger than the reaction 
plane eccentricity for the smaller Cu+Cu system. Fol- 
lowing this idea, event-by-event elliptic flow fluctuations 
have been measured and found to be consistent with the 
expected fluctuations in the initial state geometry with 
the new definition of eccentricity [46] . In this paper, 
we use this method of quantifying the initial anisotropy 
exclusively. 

Mathematically, the participant eccentricity is given as 



£2 



(3) 
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FIG. 3: Distribution of nucleons on the transverse plane for a 
y/s NN = 200 GeV Au+Au collision event with £3=0.53 from 
Glauber Monte Carlo. The nucleons in the two nuclei are 
shown in gray and black. Wounded nucleons (participants) 
are indicated as solid circles, while spectators are dotted 
circles. 



where 



and 



are the event-by-event 



(co)variances of the participant nucleon distributions 
along the transverse directions x and y (S]. If the 
coordinate system is shifted to the center of mass of the 
participating nucleons such that (x) = (y) = 0, it can be 
shown that the definition of eccentricity is equivalent to 



£2 



\J (r 2 cos(20 part )) 2 + (r 2 sin(20 part ))' 

V) 



(4) 



in this shifted frame, where r and <p pax t are the polar 
coordinate positions of participating nucleons. The 
minor axis of the ellipse defined by this region is given as 



02 



atan2 ((r 2 sm(20 part )) , (r 2 cos(2</> part ))) + n 



(5) 

Since the pressure gradients are largest along ip2, the 
collective flow is expected to be the strongest in this 
direction. The definition of v-i has conceptually changed 
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FIG. 4: Top: average elliptic flow, (172), as a function of eccentricity, £2; bottom: average triangular flow, {V3), as a function 
of triangularity, £3, in ^/s NN = 200 GeV Au+Au collisions from the AMPT model in bins of number of participating nucleons. 
Error bars indicate statistical errors. A linear fit to the data is shown. 



to refer to the second Fourier coefficient of particle 
distribution with respect to ^2 rather than the reaction 
plane 



v 2 = (cos(2(> - ^a))) 



(6) 



This change has not impacted the experimental definition 
since the directions of the reaction plane angle or tp2 are 
not a priori known. 

Drawing an analogy to eccentricity and elliptic flow, 
the initial and final triangular anisotropies can be quan- 
tified as participant triangularity, £3 , and triangular flow, 
V3, respectively: 



S3 



yj (r 2 cos(30 part )) 2 + (r 2 sin(3<?(>p ar t)) z 



(r 2 ) 



v 3 = (cos(3(> - V> 3 ))) 



(7) 
(8) 



where ^3 is the minor axis of participant triangularity 
given by 



■03 



atan2 ((r 2 sm(30 part )) , (r 2 cos(30 part ))) + tt 



3 



(9) 

It is important to note that the minor axis of trian- 
gularity is found to be uncorrelated with the reaction 
plane angle and the minor axis of eccentricity in Glauber 
Monte Carlo calculations. This implies that the average 
triangularity calculated with respect to the reaction 
plane angle or tp2 is zero. The participant triangularity 
defined in Eq. [7j however, is calculated with respect to 
■03 and is always finite. 



The distributions of eccentricity and triangularity cal- 
culated with the PHOBOS Glauber Monte Carlo imple- 
mentation [47] for Au+Au events at ^/s NN = 200 GeV are 
shown in Fig. [2} The value of triangularity is observed to 
fluctuate event-by-event and have an average magnitude 
of the same order as eccentricity. Transverse distribution 
of nucleons for a sample Monte Carlo event with a 
high value of triangularity is shown in Fig. [3j A clear 
triangular anisotropy can be seen in the region defined 
by the participating nucleons. 



IV. TRIANGULAR FLOW IN THE AMPT 
MODEL 

To assess the connection between triangularity and 
the ridge and broad away-side features in two-particle 
correlations, we study elliptic and triangular flow in 
the AMPT model. AMPT is a hybrid model which 
consists of four main components: initial conditions, 
parton cascade, string fragmentation, and A Rclativistic 
Transport Model for hadrons. The model successfully 
describes main features of the dependence of elliptic flow 
on centrality and transverse momentum |40| . Ridge and 
broad away-side features in two-particle correlations are 
also observed in the AMPT model [48l|49]. Furthermore, 
the dependence of quantitative observables such as away- 
side RMS width and away-side splitting parameter D 
on transverse momentum and reaction plane in AMPT 
reproduces the experimental results successfully, where a 
ZYAM-based elliptic flow subtraction is applied to both 
the data and the model [30] EI] . 
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FIG. 5: Dashed lines show (a) second Fourier coefficient, V2A, and (b) 
as a function of number of participating nucleons 



third Fourier coefficient, V3A, of azimuthal correlations 
200 GeV Au+Au collisions from the AMPT model. Solid 



Vpart, m yj Snn 

lines show the contribution to these coefficients from flow calculated with respect to the minor axis of (a) eccentricity and (b) 
triangularity. 



The initial conditions of AMPT are obtained from 
Heavy Ion Jet Interaction Generator (HIJING) [52] , 
HIJING uses a Glauber Model implementation that is 
similar to the PHOBOS implementation to determine 
positions of participating nucleons. It is possible to 
calculate the values of £2 , ip2 , £3 an d ipz event- by-event 
from the positions of these nucleons [see Equations [IJ [5] [T] 
and [9]. Next, we calculate the magnitudes of elliptic and 
triangular flow with respect to ip2 and ip3 respectively as 
defined in Eqs. [6] and [8] 

The average value of elliptic flow, V2, and triangular 
flow, V3, for particles in the pseudorapidity range <3 
in y/s^, = 200 GeV Au+Au collisions from AMPT are 
shown as a function of £2 and £3 in Fig. [4] for different 
ranges of number of participating nucleons. As pre- 
viously expected, the magnitude of V2 is found to be 
proportional to £2. We observe that a similar linear 
relation is also present between triangular flow and 
triangularity. 

After establishing that triangular anisotropy in initial 
collision geometry leads to a triangular anisotropy in 
particle production, we investigate the contribution of 
triangular flow to the observed ridge and broad away- 
side features in two-particle azimuthal correlations. For 
a given pseudorapidity window, the Fourier coefficients 
of two-particle azimuthal correlations, V„a, can be 
calculated in AMPT by averaging cos(nAc/>) over all 
particle pairs. Contributions from elliptic (triangular) 
flow is present in the second (third) Fourier coefficient of 
Ac/> distribution since 

J ^~2~ i l + 2v n COs(n^)} 

x {1 + 2v n cos(n(0 + Ac/)))} clcb 

= i{l + 2t£cos(nA0)}. (10) 



For a given pseudorapidity window, this contribution 
can be calculated from average elliptic (triangular) flow 
values as 

-irflow _ i £ n) 



I%M^(m){v n (vi)) Kfa))d7 7l d7 b 

I %(m)%(v2)d m d m lLl 

where n = 2 (n = 3) and the integration is over the 
pseudorapidity range of particle pairs. The average 
single-particle distribution coefficients, (v n (ji)), are used 
in this calculation to avoid contributions from non-flow 
correlations which may be present if the two-particle dis- 
tributions, v n (rji) x v n (r)2), are calculated event by event. 
The ratio (e^) / (e n ) accounts for the difference between 
(vnivt) x v n (r]2)) and (u n (?7i)) X (««(%)) expected from 
initial geometry fluctuations. 

We have calculated the magnitude of the second and 
third Fourier components of two-particle azimuthal cor- 
relations and expected contributions to these components 
from elliptic and triangular flow for particle pairs in 
y/s^ = 200 GeV Au+Au collisions from AMPT within 
the pseudorapidity range \rj\ < 3 and 2 < At] < 4. 
The results are presented in Fig. || as a function of 
number of participating nucleons. More than 80% of 
the third Fourier coefficient of azimuthal correlations 
can be accounted for by triangular flow with respect to 
the minor axis of triangularity. The difference between 
V 3 a and V3™ may be due to two different effects: 
There might be contributions from correlations other 
than triangular flow to V3A or the angle with respect 
to which the global triangular anisotropy develops might 
not be given precisely by the minor axis of triangularity 
calculated from positions of participant nucleons, i.e. 
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FIG. 6: (a) Elliptic flow, V2, and (b) triangular flow, V3, as a function of transverse momentum, pt, in bins of number of 
participating nucleons, N pa , It , for particles at mid-rapidity (|?7j < 1) in ^/s NN = 200 GeV Au+Au collisions from the AMPT 
model. Error bars indicate statistical errors. 



t>3 = ((cos(3(0 — ^3)) might be an underestimate for the 
magnitude of triangular flow. More detailed studies are 
needed to distinguish between these two effects. 

We have also studied the magnitudes of elliptic and 
triangular flow more differentially as a function of trans- 
verse momentum and number of participating nucleons 
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FIG. 7: Top: the ratio of triangular flow to elliptic flow, 
(va) I (1)2), as a function of number of participating nu- 
cleons, A^part, for particles at mid-rapidity (\rj\ < 1) in 
= 200 GeV Au+Au collisions from the AMPT model. 
Open points show different transverse momentum bins and 
the filled points show the average over all transverse momen- 
tum bins. Bottom: the ratio of different pt bins to the average 
value. Error bars indicate statistical errors. 



in the AMPT model. Figure [6] shows the results as a 
function of transverse momentum for particles at mid- 
rapidity (| ?/ 1 < 1) for different ranges of number of par- 
ticipating nucleons. The dependence of triangular flow on 
transverse momentum is observed to show similar gross 
features as elliptic flow. A more detailed comparison can 
be made by taking the ratio of triangular to elliptic flow, 
shown in Fig. [7] as a function of number of participating 
nucleons for different ranges of transverse momentum. 
The relative strength of triangular flow is observed to 
increase with centrality and transverse momentum. This 
observation is qualitatively consistent with the trends in 
experimentally measured ridge yield |25| . 



V. TRIANGULAR FLOW IN EXPERIMENTAL 
DATA 

While AMPT reproduces the expected proportionality 
of V2 and £2, the absolute magnitude of vi is underesti- 
mated compared to data and hydrodynamic calculations. 
To allow a comparison of the V3A calculations to data, 
we therefore use the ratio of the third and second Fourier 
coefficients. For data, this ratio is given by 



V 3A _ / C*(A0) cos(3A0)dA</> 
~ JC{A(j>) cos(2A0)dA<£' 



(12) 



The factors A and B in Eq. 2] cancel out in this ratio. 
Results for PHOBOS QH [25] and STAR gl] measure- 
ments are plotted as a function of number of participating 
nucleons in Figures |8(a) and |(b)| respectively. It is 
observed that V 3 a/V 2 a increases with centrality and 
with the transverse momentum of the particles. 

Also shown in Fig. [8] is the magnitude of V3A /V 2 a in 
the AMPT model with similar T], A77 and pt selections 
to the available experimental data. The calculations 
from the model show a qualitative agreement with the 
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FIG. 8: The ratio of the third to second Fourier coefficients of azimuthal correlations, V3A/V2A, as a function of number 



(a) 



of participating nucleons, 2V P art, for Au+Au collisions at ^/s NN = 200 GeV. Filled points show values derived from 
PHOBOS 1 1 91 125| and (b) STAR [41] data. Pseudorapidity and transverse momentum ranges and charge selection of particle 

Open points show results from the AMPT model for similar 
Error bars indicate statistical errors for AMPT 



pairs for different measurements are indicated on the figures, 
selection of pseudorapidity and transverse momentum to the available data, 
and combined statistical and systematic errors for the experimental data. 



data in term of the dependence of V3A/V2A on the 
pseudorapidity region, particle momenta and centrality. 
Since the V3A component of two-particle correlations 
in the model is known to be mostly due to the trian- 
gular anisotropy in the initial collision geometry, this 
observation suggests that triangular flow may play an 
important role in understanding the ridge and broad 
away-side structures in data. 

A closer look at the properties of the ridge and broad 
away-side is possible via studies of three particle correla- 
tions. Triangular flow predicts a very distinct signature 
in three particle correlation measurements. Two recent 
publications by the STAR experiment present results on 
correlations in Ac6i-Ac/>2 space for \r)\ < 1 [28] and in 
A771-A7/2 space for \A(f>\ < 0.7 [53]. In A0i-A0 2 space, 
off diagonal away-side correlations have been observed 
(e.g. first associated particle at A</>i ~ 120° and second 
associated particle at Ac/>2 « —120°) consistent with 
expectations from triangular flow. In A^-A^ space, no 
correlation structure between the two associated ridge 
particles was detected, also consistent with triangular 
flow. 



VI. SUMMARY 

We have introduced the concepts of participant trian- 
gularity and triangular flow, which quantify the triangu- 
lar anisotropy in the initial and final states of heavy-ion 
collisions. It has been shown that inclusive and triggered 



two-particle azimuthal correlations at large Arj in heavy- 
ion collisions are well described by the first three Fourier 
components. It has been demonstrated that event-by- 
event fluctuations lead to a finite triangularity value in 
Glauber Monte Carlo events and that this triangular 
anisotropy in the initial geometry leads to a triangular 
anisotropy in particle production in the AMPT model. 
The third Fourier coefficient of azimuthal correlations 
at large pseudorapidity separations have been found to 
be dominated by triangular flow in the model. We 
have studied the ratio of the third and second Fourier 
coefficients of azimuthal correlations in experimental 
data and the AMPT model as a function of centrality and 
pseudorapidity and momentum ranges of particle pairs. 
A qualitative agreement between data and model has 
been observed. This suggests that the ridge and broad 
away-side features observed in two-particle correlation 
measurements in Au+Au collisions contain a significant, 
and perhaps dominant, contribution from triangular 
flow. Our findings support previous evidence from 
measurements of the system size dependence of elliptic 
flow and elliptic flow fluctuations on the importance 
of geometric fluctuations in the initial collision region. 
Detailed studies of triangular flow can shed new light on 
the initial conditions and the collective expansion of the 
matter created in heavy-ion collisions. 
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